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Abstract Under the assumption that the photospheric quiet-Sun magnetic field 
is turbulent, the cancellation function has previously been used to estimate the 
true, resolution-independent mean, unsigned vertical flux (|B z |)truo- We show 
that the presence of network elements, noise, and seeing complicate the mea- 
surement of accurate cancellation functions and their power-law exponents k. 
Failure to exclude network elements previously led to too low estimates of both 
the cancellation exponent k and of (|-B z |)true- However, both k and (|-B z |)true are 
over-estimated due to noise in magnctograms. While no conclusive value can be 
derived with data from current instruments, our Hinode/SP results of k ;$ 0.38 
and (l-Bz |)tme ~ 270 gauss can be taken as upper bounds. 
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1. Introduction 



Turbulent fields possess self-similar power- laws, i.e., they are fractals ( Brandenburg et al., 1992| [ 
|Frisch, 1995[ ). Such power-laws (i.e, self-similarity) are found in systems where 
the underlying physical processes are the same at all scales (e.g., vortex or 
flux tube stretching suggested by IFrisch] (I1995p ). In these systems, the fields 
have the same degree of complexity (look the same) regardless of the observa- 
tional resolution. For example, incompressible magnctohydrodynamics (MHD) is 
expected ( Iroshnikov, 1964] |Kraichnan, 1965 Goldreich and Sridhar, 1995[ ) and 
found (e.o,. ILee et aU 120101 |Pietarila Graham, Mininni, and Pouquet I201ip to 
display self-similar magnetic energy spectra. Power-law magnetic energy spectra 
are also seen in compressible, stratified MHD photospheric simulations ( Brandenburg, 1995| )[ 
including those with realistic radiation and partial ionization ( |Pietarila Graham, Cameron, and Schussler, 201 
|Moll et al., 201l| . Power-laws are found in solar observations of the line-of-sight 

magnetic energy, kinetic energy, and even the granulation intensity ( Abramenko and Yurchyshyn, 2010aH 
Abramenko and Yurchyshyn, 2010b Goode et al., 20101. Power-laws are, of course ,| 
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fundamental in measurements of the fractal dimension of magnetic structures in 



both observations and simulations (e.g., Janfien, Vogler, and Kneer I2003[) . Both 



the observations and simulations above indicate that the quiet Sun magnetic field 
is likely turbulent. 

The application of turbulent self-similarity to estimate the true mean unsigned 
vertical flux from magnetograms was introduced by |Pietarila Graham, Danilovic, and Schiisslerj 
({20091 hereafter PGDS2009). Define the "net" flux (Flux;) at scale I for any 
vertical field f z to be the flux remaining after averaging over boxes A%(1) of edge 
length I: 



Fhuq=V / f z d> 



(1) 



where the boxes partition the total area: (J. >AiQ) = A. The unsigned flux of the 
field is given by 

FIuxq = f \f z \da. (2) 
J A 

The ratio of the net flux at length scale I to the true mean unsigned flux (i.e., 
flux cancellation) is called the partition function 

m Flux; (q\ 

W)=™ . 3 ) 

fTuxrj 

which follows a power-law for self-similar fields, 

x(0«r«, (4) 



where K is the cancellation exponent ( |Ott et al., 1992 ). For a completely self- 



similar field, knowledge of net flux at any scale and of the power-law scaling 
exponent implies knowledge of net flux at all scales and, hence, the unsigned 
flux. The two extreme examples are a unipolar field, for which the net flux always 
equals the unsigned flux (n = 0), and a random field. For a random field, the 
observed net flux equals the unsigned flux times the ratio of the noise correlation 
length to the diameter of the resolution element (k = 1). See the Appendix for 
a mathematical derivation of the random noise case and of the general formula, 
Equation (fl"5"j) . to determine the unsigned flux for a completely self-similar field 
given k and the correlation length below which the field becomes smooth. 

Under the assumption that the quiet-Sun magnetic field is turbulent, the 
theoretical framework in the Appendix applies. This might suggest that the 
unsigned quiet-Sun flux could be deduced from Zeeman-based instruments. Such 
studies have previously been made (PGDS2009: IStenflol l2TJTTj> . PGDS2009 ob- 
served a power-law in x(l) (which they dubbed the "cancellation function" when 
applied to a magnctogram) down to « 200 km with k = 0.26 in Hinode Spectro- 
Polarimcter (SP: ITsuneta et aTll2008p data. Note that power-laws for \(P) have 
previously been seen in simulations of the electric currents ( |Sorriso-Valvo et al, 2002| | 



Pietarila Graham, Mininni, and Pouquet, 2005 ) and in solar observations of current -| 
helicity ( |Sorriso-Valvo et al, 2003[|Abramenko,"2003{|Sorriso-Valvo et al, 2004| 



|Abramenko and Y urchyshyn, 2010a). PGDS2009 used Equation (fl"5"j) to estimate 
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the net flux at a scale of 800 m to be sa 50 G. Turbulent power-laws (e.g., energy 
spectra and extend only down to the dissipative range, but the magnetic 
diffusion scale is expected to be significantly smaller than 800 m. Thus, the 
extrapolation was taken as a lower bound, Fluxo > 50 G (gauss). This result 
was found to be in agreement with a separate extrapolation from radiative MHD 
simulations. 

Other estimates of k can be made. ILites et~a l. (2008) found a 30% increase 
in flux changing from a resolution of 1" to 0.33" which indicates k — 0.24. 
Sola nki et al.\ (|2010p found a factor of 3.7 more flux on doubling of resolution, 
indicating k — 1.9. Using many instruments with different spatial resolutions, 
San chez Almeida and Martine z Gonzalezl (|201ip fit Flux/ oc i.e., k=1. The 
data set used by PGDS2009 was recalibrated bv lStenflol (|201ip leading to more 
pronounced high field strength tails in the flux probability distribution. The 
power law from the recalibrated magnetogram yielded a value of n that is half 
of the value derived by PGDS2009. In addition to the SP magnetogram [Stenflol 
(2011) computed cancellation functions from a set of Hinode Narrow-band Filter 
Imager (NFI) Na I D\ magnetograms. The resulting k was found to be 0.127, 

1. e., very similar to the value for the recalibrated SP magnetogram. In order to 
agree with Hanle-based magnetic field measurements, this cancellation function 
needs to be extrapolated down to a spatial scale of 10-100 m. 

Another application for using the fractal nature of magnetic fields is in flare 
predictions. Previous works (|Sorriso-Valvo et aU l2~003. 200jl |Abramenko and Yurchyshynj [ 
I2010ap have shown it to be promising. However, as in the case of cancellation 
functions, the data-sample sizes used to study the suitability of the method have 
been limited. Recently |GeorgoinTs| (|2012p has shown how previous analyses may 
have been misleading: A statistical comparison of flaring and non-flaring active 
regions found that the fractal properties of the linc-of-sight magnetic field cannot 
distinguish flaring and non-flaring active regions. Many of the previous works, 
however, used current helicity, not magnetic flux, for the analysis and therefore 
the conclusions of Georgoulis (2012) may not apply to them. 

Prompted by the results of PGDS2009,lStenf]o|([lIEI]), and |GeorgouHsl flUED) , 
we suggest that measurements of the cancellation function do not confront the 
reality of the observational data. Such confrontation is the aim of this paper. 
We analyze the statistics of quiet Sun cancellation functions using magnetograms 
from three different instruments. We test if the cancellation exponent k is robust 
within and among the instruments to point out observational artifacts in the 
measurement of n and to determine if any conclusive values or bounds can be 
made with existing data. 

2. Data and Methods 

We use altogether «800 magnetograms from the Helioseismic and Magnetic 
Imager (HMI: fScherrer et aU I2012"f on the Solar Dynamics Observatory (SDO) 
satellite, Michelson Doppler Imager (MDI; [Scherrer et aUU995)l on the Solar and 
Hcliosphcric Observatory (SOHO) sarellite, and Hinode SP. The largest data set, 
700 HMI magnetograms, allows us to characterize the variation within a single 
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instrument. Data from MDI and SP allow for inter- instrument comparisons and 
for identification of possible instrumental artifacts as well as for determination 
of the effects of spatial resolution and magnctogram noise. 

The HMI data consist of daily (from end of March 2010 to mid- June 2011) 4 s 
integration magnetograms taken around 12:00 UT. We analyze a 501 x 501 pixel 
(w 253 "x253 ") area around the disk center. For MDI, we use 301 x 301 pixel 
(w 182 "xl82 ") sub-regions, excluding active regions, in high resolution magne- 
tograms (level 1.8) near the solar disk center. Note that a calibration coefficient 
(BSCALE=2.81) is applied to the MDI magnetograms to convert from counts to 
gauss leading to the magnetograms being quantized in units of BSCALE. 

The SP magnetograms are the longitudinal magnetic field measurements in 
Level ID data with exposure times varying between 1.6 and 12 s. Both «0.16" 
and «0.32" pixel size magnetograms are included in the data set. We choose ob- 
servations of quiet Sun taken near the disk center. The size of SP magnetograms 
ranges from w 400 x 400 to 1000 x 1000 pixels (« 64 "x64 "-160 "xl60 ")• 

Cancellation functions [x(0] are computed for each magnetogram using the 
Monte Carlo technique of ICadavid et al.\ ([1994). Linear fits to x(Z) in log-log 
space are made for three different spatial scales: small scale <2" (exponent Ks), 
intermediate scale 2-5" (k\) and large scale 2-9 " (kl)- While the physical 
scales are the same for all instruments, the number of pixels for fitting the 
exponent varies due to the instruments' different pixel sizes (HMI «0.5", MDI 
«0.6", SP «0.16"or «0.32"). The spatial scales were chosen to be sensitive to 
the scales most affected by noise (ks) and scales not dominated by noise, but 
still below the spatial scale of the network (ki and kl). The overlap of Ki and 
kl ensures that enough data points (14 for HMI, 11 for MDI and 24 or 47, 
depending on pixel size, for SP) are included in fitting kl and that the results 
are not strongly sensitive to the upper cut-off (5" or 9") of the fitted spatial 
scales. 

Additionally, the mean unsigned flux (\B Z \) and flux imbalance (B z ) / (\B Z \) 
are computed for each magnctogram. 

3. Results 

3.1. Exponents and Effects of Network 

The results of the analysis on the HMI cancellation function are summarized in 
Figure [T]and Table [TJ Compared to large spatial scales (>10", for which we do 
not fit a power- law), the exponents at the smaller scales are fairly similar: kt and 
KLare nearly identical, which indicates that we are fitting a genuine power-law 
at these scales ( 2"to 9"). \(P) turns up at the small scales as is demonstrated 
by ks having ps30% higher values than k,\ and kl- The increase is partially 
due to contributions from noise becoming more dominant (see Section 13.31 and 
note that K noiso = 1: ISorriso-Valvo et al.\ (|2002|) and the Appendix). Note that 
there is both scatter (standard deviation up to 6% of the mean) and systematic 
differences between the exponents of power-law fit to different magnetograms. 

All exponents decrease as a function of (|B 2 |), (\B Z \)= Flux; for I = 1 pixel. 
The more flux there is in the magnetogram, the smaller each k is. A linear fit 
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Figure 1. HMI cancellation functions. Left: Cancellation functions for selected magnetograms 
(in various colors) normalized to unity at 9". Vertical and horizontal lines mark ranges used for 
fitting the exponents (Kg in green dotted line and crosses, n\ in red dashed line and diamonds, 
and kl in black dash-dotted line and circles). The same colors and symbols are used in all 
plots to differentiate between the spatial scales. Middle: Cancellation exponents as a function 
of (|S Z |). Error bars show lufit- Right: Exponents as a function of flux imbalance. 



Of function of log(|£? 2 |) shows that the decrease is strongest for kl and 

weakest for /eg. No correlation is found for k and the global (full-disk) unsigned 
flux, which reflects the overall activity present on the solar disk at a given time: 
The cancellation functions reflect the local magnetic environment, whose fractal 
properties do not appear to be influenced by the global field, /e does not depend 
on the flux imbalance at small imbalances (below mlO%). For imbalances greater 
than «10% k decreases with increasing imbalance on all spatial scales. 

The exponents show no strong temporal evolution or indications of the ex- 
ponent changing as solar activity increases (Figure [2]). Since the activity in the 
rising phase of the solar cycle emerges at high latitudes, it is not surprising that 
no significant change in the exponent is seen near disk center. 

The results for MDI are shown in Figure [3] and Table [T] They are qualita- 
tively similar to HMI, confirming that the trends ((\B Z \), flux imbalance) are of 
solar, not instrumental, origin. The exponents are larger and have more scatter 
(standard deviation is 9% of the mean). Also the dependence of k on (\B Z \) is 
stronger. 

The SP magnetograms (Figure 2] and Table [1} are a mixture of data with 
different pixel sizes and exposure times. A comparison of Ks for magnetograms 
with exposure times less than two seconds and ones with above, shows that /eg 
is larger for shorter exposure (and thus noisier) magnetograms. The difference 
is smaller for K\ and kl- This is to be expected, if the difference is due to 
noise, which is more dominant at small scales (see Section [575)1 . The measured 
exponents are similar to k found by PGDS2009 who measured /t=0.26 from 
a single magnctogram. Regardless of exposure time, all the SP exponents are 
smaller than MDI and HMI exponents. No dependence of k on (\B Z \) or flux 
imbalance is seen. The sample, however, is too small to establish this. Note 
that for the magnetograms with the longest exposure times a mixture of spatial 
and temporal cancellation may occur. The time to raster 2" varies strongly 
depending on the exposure time and step size: For exposure times of 1.6 and 4.8 
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Figure 2. HMI cancellation exponents as a function of time. Kg is shown in green crosses, K\ 
in red diamonds, and kl in black circles. (|£J z |)is shown in dark blue * (y-axis on right). 
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Figure 3. MDI cancellation functions. Left: Cancellation functions for selected magnetograms 
(in various colors) normalized to unity at 9". Vertical and horizontal lines mark ranges used for 
fitting the exponents (Kg in green dotted line and crosses, kj in red dashed line and diamonds, 
and kl in black dash-dotted line and circles). The same colors and symbols are used in all 
plots to differentiate between the scales. Middle: Cancellation exponents as a function of (|i? z |). 
Error bars show lo"fl t . Right: Exponents as a function of flux imbalance. 



s the rastering times are well below the granulation turnover time scale. For the 
exposure time of 8 s and above the rastering time is over 100 s. 

The dependence of MDI and HMI cancellation exponents on (|-B Z |) and flux 
imbalance suggest that the presence of magnetic network elements in the magne- 
tograms influences the cancellation exponents. (Magnetic network flux is stronger 
than the intra-network, the network is mostly unipolar, and occurs at scales 
of a few tens of arcseconds.) To model the effect on x(0 of a large (active- 
region-remnant or network-like) unipolar flux region, our observation region A 
is decomposed into two sub-regions: B (containing intra-network quiet-Sun in 
which a scaling %(/) would be found) and the unipolar flux concentration C 
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Table 1. Power-law fits to x(0 for HMI, MDI, and SP magnetograms. 





Mean 


Standard deviation 




i 


HMI 2 










«S 


0.53 


0.017 


0.74 (0.01) 


-0.22 (0.01) 


«i 


0.36 


0.020 


0.61 (0.01) 


-0.26 (0.01) 




0.35 


0.022 


0.66 (0.01) 


-0.34 (0.01) 


MDI 2 










res 


0.61 


0.053 


1.62 (0.06) 


-0.64 (0.03) 


«i 


0.45 


0.033 


1.47 (0.07) 


-0.65 (0.04) 




0.42 


0.036 


1.51 (0.08) 


-0.70 (0.05) 


SP 










^exp ^2 S 












0.30 


0.020 








0.28 


0.020 








0.29 


0.014 






^exp ^2 S 










«s 


0.24 


0.0063 








0.26 


0.014 








0.28 


0.032 







1 K = c+/*log((|B z I)) CTfl t is given in parentheses. 

2 Only magnetograms with flux imbalances below 5% are included. 
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Figure 4. SP cancellation functions. Left: Cancellation functions for selected magnetograms 
(in various colors). The cancellation functions with pixel size 0.32"are normalized to unity 
at 18"and pixel size 0.16"at 9". Vertical and horizontal lines mark ranges used for fitting 
the exponents (reg in green dotted line, k\ in red dashed line, and kl in black dash-dotted 
line). The same colors are used in all plots to differentiate between the spatial scales. Middle: 
Cancellation exponents as a function of (|B^|). Error bars show l<Tf} t . Circle symbols arc rasters 
with exposure times less than 2 s, diamonds are exposure times between 2 and 5 s, and squares 
are exposure times longer than 5 s. Right: Exponents as a function of flux imbalance. 
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(over which J c \B z \da = | j c B z da\, i.e., there is no cancellation). Our goal is 
to find the cancellation function [x'(0] over the entire region A = B U C. We 
consider only scales I <C the size of the unipolar flux concentration. So, a given 
square Ai (I) is either in B or C, 



X'd) = 



E 



E 



Ai(i)ec 



Jam BAa 



Using EU(i)GC 



J B \B z \da + J c \B z \da 
= J c \B z \da and defining 



P = J c \B z \da/ Jg \B z \da, the ratio of mean unsigned flux in the unipolar flux 
concentration compared to the rest of the region, we find 



x'(l) = 



x(i) + P 
l + p 



(0) 



Assuming that the intra-network quiet-Sun flux is completely balanced, flux 
imbalance (FI) is given by 



Fl = 



J A B z da 



! A \B z \da 1 + P 



(7) 



We take x(7) from the HMI magnetogram with the smallest flux imbalance 
and employ Equations ([5]) and ([7]) to plot synthetic x'(0 versus flux imbalance 
in Figure [5] The synthetic exponents as a function of flux imbalance mimic 
the MDI and HMI observations: Small increases in the imbalance do not alter 
the exponent, while imbalances above «10% lead to a decreasing k. The model 
explains the observed dependence of k on flux imbalance as a consequence of 
unipolar network magnetic fields altering the cancellation function. Note that 
the dependence of the exponents on (\B Z \) applies also to magnetograms with 
very small imbalances. Small imbalances imply that the positive and negative 
network concentrations balance out each other, not that there is only a small 
amount of network flux present. 



3.2. Lntra-Network Exponents 



Since the cancellation exponents are affected by the presence of network ele- 
ments, a combination of smoothing and thresholding can be applied to mask out 
all the network pixels in the magnetograms prior to computing the cancellation 
functions. An example of a mask applied to an HMI magnetogram is shown in 
Figure |U To test how the masking affects the measured exponents, we apply a 
network mask to a synthetic magnetogram consisting of pure noise (K no isc = I ) 
and compare the exponents from the masked and unmasked magnetograms: The 
effect on the exponent is smaller than the error in the fit. 

Masking out the network pixels in the HMI magnetograms to measure the 
"true" intra-network cancellation exponents (Figure [7] and Table [5]) leads to 
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Figure 5. Effect of flux imbalance. Synthetic cancellation exponents for intermediate scale 
Kj as a function of flux imbalance. Error bars are l<7fi t . Measured HMI kj are shown as red 
crosses. 




Figure 6. Masking of network. Left: Original HMI magnetogram with network elements. 
Right: Same magnetogram with network elements masked out. The color scales saturate at 
±25 G. The horizontal line in the bottom left corner shows the spatial scale of 50". 
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Figure 7. HMI cancellation functions for intra-network (network masked out) magnetograms. 
Left: Cancellation functions for selected magnetograms (in various colors) normalized to unity 
at 9". Vertical and horizontal lines mark ranges used for fitting the exponents (res m green 
dotted line and crosses, K\ in red dashed line and diamonds, and «l m black dash-dotted line 
and circles). The same colors and symbols are used in all plots to differentiate between the 
spatial scales. Middle: Cancellation exponents as a function of Solid black line shows 

linear fit, slope 0.3, of k vs. \og(\B z \). Error bars show l<Tfi t . Right: Exponents as a function 
of flux imbalance. The systematic dependence is seen to be effectively removed. 

significantly higher values for kl and K\. k still depends on (\B Z \), but in 
the opposite sense than in the unmasked magnetograms: In the intra-network 
magnetograms k increases with increasing (|-B 2 |). This demonstrates that the 
decrease of n with increasing ( | B z | ) in the unmasked magnetograms was due to 
the presence of the network elements. A linear fit of intra-network magnetogram 
k vs. \og{(\B z \)) gives a slope of «0.3. This increase of k with (|-B Z |) can be 
understood from Equation Q . An increase in flux imbalance means more pixels 
that arc never canceled out at any scale, a flattening of the cancellation function, 
and a decrease in k. A magnetogram with very weak signal can be out of balance 
from a small number of moderately strong pixels of the same sign. In fact, we 
see in Figure [5] that nearly all magnetograms with over 1% flux imbalance, have 
(|-B 2 |}< 5.5 G. The weakest magnetograms have greater flux imbalances and, 
consequently, smaller k. 

The network masking of MDI magnetograms was not entirely successful and 
effects of network are still visible in the (|-B Z |) and flux imbalance plots. Com- 
pared to HMI, the effect of masking the magnetograms is small in MDI (Figure 
[5] and Tabled]). The modest change may also be partially due to the network 
vs. intra-network contrast being smaller in MDI than HMI. The intra-network 
SP exponents (Figure fTTJI and Table (2) are increased on all scales. The change is 
largest for longer exposures and large spatial scales. Since the longer exposure 
times have less noise and their k values increase the most after removal of the 
network, the increase in k cannot be due to noise alone. 

3.3. Effect of Noise and Seeing 

An examination of Tables [T] and [5] suggests that the cancellation exponent is 
larger for noisier instruments and shorter exposures, i.e., for noisier magne- 
tograms. We now demonstrate that this is the case by artificially adding noise 
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Figure 8. Flux imbalance in masked HMI magnctograms vs. (|B Z |). 
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Figure 9. MDI intra-network (network masked out) cancellation functions. Left: Cancellation 
functions for selected magnctograms (in various colors) normalized to unity at 9". Vertical and 
horizontal lines mark ranges used for fitting the exponents (Kg in green dotted line and crosses, 
Ki in red dashed line and diamonds, and kl in black dash-dotted line and circles). The same 
colors and symbols are used in all plots to differentiate between the scales. Middle: Cancellation 
exponents as a function of (|B Z |). Error bars show l<7f} t . Right: Exponents as a function of 
flux imbalance. 



to our magnetograms. The first panel in Figure [TT] shows how an HMI can- 
cellation function changes as random-distributed noise is incrementally added 
to the magnetogram prior to computing %(/). Adding noise changes first the 
smallest scales. As the amount of noise increases the scales affected by noise 
also increase. For HMI adding I G (random numbers with a mean of zero and 
a standard deviation of one) of noise does not change the cancellation function 
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Figure 10. SP intra-network (network masked out) cancellation functions. Left: Cancellation 
functions for selected magnctograms (in various colors). The cancellation functions with pixel 
size 0.32"are normalized to unity at 18"and pixel size 0.16" at 9". Vertical and horizontal lines 
mark ranges used for fitting the exponents (kq in green dotted line, kt in red dashed line, and 
kl in black dash-dotted line) . The same colors are used in all plots to differentiate between the 
spatial scales. Middle: Cancellation exponents as a function of (|£J Z |). Error bars show lcrflf 
Circle symbols are rasters with exposure times less than 2 s, diamonds are exposure times 
between 2 and 5 s, and squares arc exposure times longer than 5 s. Most of the scatter is due 
to differing exposure times. Right: Exponents as a function of flux imbalance. 



Table 2. Power-law fits for intra-net- 
work (network masked out) cancellation 
functions for HMI, MDI, and SP magne- 
tograms. 





Mean 


Standard deviation 


HMI 










0.63 





024 
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MDI 1 








«S 


0.59 
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020 


><i 
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«S 
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><l 


0.38 





014 


KL 


0.48 
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1 On\y magnctograms with no network remnants are included. 
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Figure 11. Noise and x(0- Left: HMI intra-network cancellation functions with increasing 
artificial noise (see figure legend). Middle: Histograms of intra-network magnetic flux. Right: 
Small-scale exponents as a function of added noise. For the SP magnctogram we use a raster 
taken with a 1.6 s exposure time. 



noticeably and only a small difference is seen for 3 G. In contrast, adding noise at 
the 6 and 15 G levels, respectively, significantly alters on increasing spatial 
scales. That no change is visible for 3 G, gives us a measure of noise in the HMI 
magnctogram: between 3 and 6G. 

The remaining panels in Figure [11] show how the small scale cancellation 
exponents for the different instruments react to increased noise. (Noise was 
added in units of counts (1 count = 2.81 G) to the MDI magnetograms.) All 
of the Kg, as a function of noise curves, have the same general shape (rightmost 
panel in Figure [TTj) : They start with a plateau and after an instrument-specific 
threshold (magnctogram noise level) is reached, res begins to increase linearly as 
a function of added noise. The nominal noise values (upper limits) of the different 
instruments can be defined as the full width at half maximum of Gaussian fits to 
the magnetic flux histograms. The widths are 31.6 G for MDI, 21.1 G for HMI, 
and 13.5 G for SP. In the linear increase regime the increase in k per added 
noise is 0.01 per G for MDI and 0.02 per G for SP and HMI. The length of the 
plateau and steepness of the increase are related to the inherent noise in the 
data: The less instrumental noise there is, the more sensitive the cancellation 
function is to added noise. SP, which has the lowest nominal noise level, has the 
shortest plateau and steepest increase, while MDI, the noisiest instrument, is 
less sensitive to noise. Consistent with the noise experiment, SP has the smallest 
measured exponents and MDI the largest, suggesting that MDI magnetograms 
are more strongly affected by noise. 

Figure [T^] shows how spatial smoothing changes the cancellation function. The 
effect of smoothing is opposite to that of increasing the noise (recall that k = 
for a smooth field and k = 1 for noise). The exponents decrease with increased 
spatial smoothing. This demonstrates how seeing conditions and spatial over- 
sampling can reduce the measured exponents. 
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Figure 12. Effects of spatial smoothing. Left: HMI cancellation functions computed with 
increased spatial smoothing (Gaussian smoothing with full width at half maximum (FWHM) 
increasing in increments of 0.5 pixels). Right: Exponents as a function of Gaussian FWHM 
(«S m green, k\ in red and kl in black). 



4. Conclusions 

Significant artifacts in measuring the cancellation function from current obser- 
vations have been identified. The cancellation exponent is sensitive to noise level 
differences between instruments: Noise increases both k and the scatter between 
individual measurements of n. Not surprisingly, k decreases by seeing and spatial 
over-sampling. We also find that k is sensitive within a single instrument's ob- 
servations to magnetic flux imbalances and mean unsigned flux. A simple model 
suggests that this dependence is the effect of network elements: They reduce k. 
This simple model also explains why the value of k calculated in PGDS2009 and 
Ste nflol (|201f p differ: The network is enhanced in the latter study, decreasing k. 

Of all the observational artifacts, the effect of the network can be removed 
by masking network elements out of magnetograms. We then find values of 
ki « 0.58 for HMI and kj « 0.38 for SP. (For MDI the removal of network 
was not fully successful and the dependence of x(l) on the presence of network 
elements was not entirely removed.) The difference in k between HMI and SP is 
due to HMI magnetograms being demonstrably more noisy: HMFs noise level is 
« 4 G and SP's is ps 2 G, at least for determinations of the cancellation exponent 
(see Figure HI]). Since noise first affects the smallest scales of the cancellation 
function, high resolution (diffraction limit, pixel size, and seeing) observations 
with minimal noise are needed to measure the exponent on scales below a couple 
of arcseconds. Therefore, no conclusive value of the unsigned flux can be given 
by cancellation function analysis with current instruments. 

An upper bound of the unsigned flux is the best that can be done with 
the available data. Since SP has the lowest noise level, we can safely conclude 
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k ;$ 0.38. The quiet-Sun magnetic field is not a completely self-similar field: 
It will become smooth over an order-of-magnitude in scales near the magnetic 
diffusion scale (below which there is no further cancellation). However, using 
k ^ 0.38 and that the power-law ends at Iq > 80 m (the magnetic diffusion 
scale estimate from PGDS2009) in Equation (fT5|) . we have an upper bound 
on the true mean unsigned flux. We find that (|-B z |)true = Fluxo 268 G 
(standard deviation 44 G) from SP data with exposure times longer than 5 s. 
(Note that if the power-law for some reason ends at even larger scales than the 
« 800 m beginning of the diffusive range, this bound still holds.) The noisier 
instruments, MDI and HMI, give significantly higher values. Both the noise 
level and sensitivity of the instruments affect the estimates. SP being the least 
noisy of the instruments can be considered as an upper bound for the flux. 
|Manso Sainz, Lan di Dcgl Tnnocenti, and Trujillo Bueno| (|2006[) found a similar 
upper bound, (|-B z |)tmc ~ 200 — 300 G from the scattering line polarization of 
Ce II. 

We found no correlation between k and global mean unsigned flux: The small 
scale (<10") field distribution/fractal geometry is not affected by the global 
magnetic configuration, at least not during a rising phase of the cycle sampled 
by HMI. A study of HMI cancellation functions from the end of the previous 
solar minimum past the next solar maximum may show how flux from decaying 
active regions affects the flux distribution in the quiet Sun, and possibly give 
indications of the relative importance of possible quiet-Sun small-scale dynamo 
action in different phases of the solar activity cycle. 

It should be emphasized that the present analysis and findings are for quiet 
Sun magnetic fields which are known to be turbulent. Issues identified in the anal- 
ysis most likely do not affect as strongly measurements of active region fractal 
dimensions. The significance of noise in the cancellation statistics may extend 
to the studies of fractal dimensions and flaring probability (e.g., |Georgoulis| 
I2012p . Based on the current analysis, using less-noisy magnetograms, such as 
HMI, may show the fractal dimension of the magnetic flux to still have some 
predictive capability for flaring active regions. To establish the usability of fractal 
analysis for flare forecasting, however, a large sample of flaring and non-flaring 
active region magnetograms with little noise is needed. A statistical study of 
flaring and non-flaring HMI vector magnetograms could also be used to address 
the differences of using line of sight flux and current helicity for predictions. 
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Appendix 

For a completely self-similar field, knowledge of net flux at any scale and of 
the power-law scaling exponent implies knowledge of net flux at all scales and, 
hence, the unsigned flux. This is most easily seen in the example of a vertical 
random field [f z ]. For a random field, the integral over an area is proportional 
to the square root of the area (random walk) , 

f z da cxZ, (8) 

while the integral over a strictly positive field is proportional to its area, 

|/ z |dacxZ 2 . (9) 



L 



Therefore, 

x(0«j, (io) 

and we identify k = 1 for pure noise. For scales smaller than the correlation 
length of the noise, the net flux is equal to the unsigned flux, 

Flux/ = Flux V/ < l , (11) 

and for larger scales, 

Flux/ = — Vl>l - (12) 

The constant C is determined by equating Equations (fTTj) and (fT2")l at the 
correlation length, 1 = 1$: 

Flux* l o w , . , riq -. 

Flu^ = 7 V/ " V (13) 

This power-law dependence, Equation. (|f 3I) . is shown in Figure [T5l This means 
that the unsigned flux can be exactly calculated given the net flux at any scale 
and the scale at the end of the power-law, Iq, 

Flux = Flux/-^-, (14) 
*o 

for a random field. For any other purely self-similar field with cancellation 
exponent k, 

Flux/ = Flux L (y) K VZ, L > l ■ (15) 

For the extreme case of a unipolar field, k = and Flux/ = Fluxo for all scales 
I. In general, < k < 1. 



SOLA: ms_spl.tex; 1 October 2012; 0:37; p. 16 



Cancellation Functions of Quiet Sun Magnetic Flux 



17 




1 10 100 

1/IJ0 

Figure 13. Ratio of net flux to unsigned flux versus scale (I) for noise with correlation 
length Zq- The ratio is given by -f- (dashed line). 
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